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Abstract 

We express the leading electromagnetic corrections in K tttt as in- 
tegrals over the virtual photon squared-momentum Q^. The high 
behavior is obtained via the operator product expansion. The low 
behavior is calculated using chiral perturbation theory. We model 
the intermediate region using resonance contributions in order to 
enforce the matching of these two regimes. Our results confirm our 
previous estimates that the electromagnetic corrections provide a rea- 
sonably small shift in the A J = 3/2 amplitude. 
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Figure 1: Leading electromagnetic correction to K ^ vrvr. 

1 Introduction 

In a previous publication [Q] , we calculated the leading electromagnetic cor- 
rections to — > TTvr nonleptonic decays within chiral perturbation the- 
ory (ChPT).|^ The only hadronic degrees of freedom were the pseudoscalar 
mesons. Loop integrals were analyzed in terms of dimensional regulariza- 
tion, and counterterm amplitudes were introduced to cancel all divergences. 
The finite counterterms parameterize the short distance effects of heavy de- 
grees of freedom. Our ChPT analysis yielded effects which were estimated 
to be at the several per cent level. Unfortunately, due to the presence of 
many unknown finite counterterms the results were accompanied by error 
bars as large as the signals. 

1.1 The Method of Dispersive Matching 

In this paper we extend the previous calculation to higher energies by using 
a 'dispersive matching' approach. Here, active degrees of freedom include 
not only the ground state pseudoscalar mesons but also the spin-zero and 
spin-one low-lying meson resonances. Amplitudes are expressed as integrals 
over the virtual photon euclidean squared-momentum . Within chiral per- 
turbation theory, the integral is regulated dimensionally, and unknown 
constants are introduced to parameterize the contributions from intermedi- 
ate and high energy. In contrast, the dispersive matching approach is an 
attempt to construct an intermediate energy contribution that sucessfully 
interpolates between the low and high energy regions. This allows the full 
integral to be calculated. 

The K — > TTTT amplitude with EM interactions present is given generally 

'^By 'leading' is meant the component which arises from electromagnetic corrections to 
the (large) AJ = 1/2 amplitude (c/Fig. |l]). This approach will be followed here. 
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Figure 2: High-Q^ electroweak dynamics of quarks. 
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where D^^^[q) is the photon propagator and W^^{q,p) describes the scat- 
tering 'jK{p) 77r7r. Rotation to euchdean momentum space followed by 
evaluation of the angular integral yields the equivalent form 



To determine the K — > vrvr amplitudes to order e^p"^ in chiral power counting 
requires knowing Wf^^ {q,p) (or equivalently W (Q^)) at order p'^ and for all 
values of Q^. We have rigorous information on W^j^iy il^p) only in the two 
asymptotic regimes of (i) low Q^, where ChPT provides the appropriate 
couplings, and (ii) high Q^, where the quark degrees of freedom couple to 
the photon according to the standard electroweak theory. Our goal will be 
to match these regions. 

Consider the process of building W (Q^) from the low-Q'^ end. In prin- 
ciple we can use ChPT to generate W^i, {q,p) to order p^Q^". At the lowest 
energies, the dominant contribution is from the ground state mesons. At 
intermediate energies resonance degrees of freedom become active, and we 
can use effective lagrangians to describe their interactions. Such resonance 
contributions serve to soften the polynomial behavior as increases. [Q] 
Eventually, the low/intermediate description is matched to the high-Q^ 
effects of Fig. ^. 

Neither of the short distance contributions depicted in Fig. 2 plays a 
dynamical role in the radiative problem but for different reasons. The pro- 
cess of Fig. |2|(a) leads merely to an overall shift in the strength of the weak 
interaction but does not give rise to mixing between the isospin amplitudes. 
The electroweak penguin operators of Fig. ^(b) do contribute to K — > vrvr 





(2) 



2 



decay, but are found to be quite small [Q] and so are neglected in the work 
reported here. Physics of the low-to-intermediate region is therefore the 
dominating influence in our calculation. 

In Section 2, we define the various interaction lagrangians which are 
needed in the course of the calculation. We present a detailed account of the 
calculational program in Section 3, from its content through to the results 
and some phenomenological implications. We pay particular attention to the 
uncertainties inherent in our calculation, and attempt to provide realistic 
error estimates. Final remarks appear in Section 4. 

2 Effective Lagrangians 

Our starting point will be a tree level calculation of the {W-^^ {Q'^)} including 
as intermediate states the ground state mesons and the low lying resonances. 
Their interactions are dictated by the lowest order chiral lagrangians (of 
order p^). Specifically, in the resonance sector we include the vector (V), 
axialvector (A), scalar (S) and pseudoscalar (P) octets and the scalar (SI) 
and pseudoscalar (PI) singlets. 

2.1 Ground State Mesons 

The lASI = 1 octet lagrangian which governs the spinless ground state 
mesons begins at chiral order p^, 

Cf = 58 Tr (X^D^UD^'U^) , (3) 

with ^(g — 6.7 • 10^^ and U = exp [iX ■ <1>). The corresponding AS* = 
strong/electromagnetic lagrangian is 

42 = ^ Tr (d.UD^^U^) +^Tt (xC/t + Ux^) , (4) 

where x = 2i?o diag{mu,md,ms). and D^C/ = 3^C/ + ie[Q,U]A^, with 
being the photon field. 
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2.2 Spin-One Resonances 

The spin-one vector and axialvector resonances which enter our calculation 
are represented respectively by the field matrices Rfj_y = V^u , A^i, , 



and 



A 



K*- 



iyV2 + h/VQ 



yV2 + uJ8/VQ 



K*+ 



Kt 
-2/1/^6 



The normalization of Rnu is given by 



(5) 



(6) 



Mr 



(7) 



Analogous to interactions among the spinless ground-state mesons, inter- 
actions of the resonances are likewise given in terms of effective lagrangians. | 
For AS* = vertices we have 



C. 



iS = ^ Tr ( ) +^^Tr {V.^u^u^ + ^ TV {A.^f^ 



where 



U = uu 



lU 



The couplings Fy , Gy , Fa have the numerical values 0] 



(8) 



(9) 



F^ / 
Fy ~ 0.154 GeV , ~ -f , Fa ~ (F, 

Fy V 



2 p2 



1/2 



(10) 



Although the effective lagrangian used to describe \^S\ = 1 interactions of 
the resonances is given most generally by ||5| 



10 

E 

k=l 



(11) 
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only four of the ten possible operators are relevant to our K vrvr analysis, 

i^p) = Tr (A[i?^,,/r]+) , i^f ) = Tr {A[R,,jn+) , 

=iTT {A[R^„[u^^,u']U) , Kf'^=iTT{Au^R^,u,) , 

where A = wX^u. Use of the {Kj^ } introduces eight couplings {gj^ } 
{R = V,A and k = 1,2, 5, 6) into the calculation. It is convenient to convert 
these to dimensionless quantities, 

.r=^5r, 4*'=^#'. 1,2.5,6) . (13) 

2.3 Spin-zero Resonances 

Finally, we list effective lagrangians for the spinless resonances, including 
the octet scalars S, the singlet scalar and their pseudoscalar analogs P 
and Pi. We begin with the strong lagrangians, 

= CdTr{Su^un+CmTr{Sx+) 

+ CdSi Tr (u^u/^) + CmSi Tr (x+) , 

jC^r''"^ = idmTr{Px-)+idmPiTr{x-) , (14) 

where x± = ux^u ± tx^xu^. The weak lagrangian for the octet spinless 
resonances is 

4r*^ = E 9hKi + E ghKf , (15) 
i=l i=l 

where 

=Tr{A [S, x+]+) , Ki=Tr{SA)-Tr (x+) , 

i^f = Tr (A [S, X-]+) , K!=Tt (A [5, n^7i^]+) , (16) 

Ki = Tr {Au^') ■ Tr (u^S) , = Tr (AS) - Tr (ui^uf) , 

and 

K[ = i Tr (A [P, x-]+) , i^2^ = i TV (A P) ■ Tr (x-) , . . 

Kl' = z Tr (A [x+,P]) , K[ = iTriA [P, u^.u'']) . ^ 

The weak lagrangian for the singlets is given by 

^(si^glet) ^ -l^^P + ^ , (18) 

1=1 

with 

i^f = 5i TV (A x+) , i^2^ = 5i Tr (A u^.u'^) , K[ = iPi Tr (A x-) • 

(19) 
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(a) (b) (c) 

Figure 3: Contributions to Wj^'^: (a) Born, (b) resonance, (c) loop. 

3 Details of the calculation 

We have at hand the tools to construct reliable expressions at low and in- 
termediate for the {VFj(Q^)} functions of Eq. @. The two major com- 
ponents will be: 

1. Tree diagrams (Figs. ^a),(b)) involving exchange of the ground state 
pseudoscalar mesons (Born terms) and of resonances: 

Within chiral perturbation theory, the vertices in the tree dia- 
grams are described by point-like couplings at leading order. However, 
in QCD we know that the couplings fall off at higher Q^. In order to 
incorporate this feature, we model form factor corrections to the Born 
terms with vector resonance contributions. The set of Born diagrams 
(together with insertions of meson form factors) is free of unknown 
parameters. 

The remaining vector and axialvector resonance contributions de- 
pend on eight unknown weak couplings. Various phenomenological 
inputs can be used to fix them, but some remain unconstrained. In 
principle this part of the amplitude requires matching to the penguin 
short distance contribution. The requirement that matching occurs 
successfully affords a way to further constrain the unknowns. This is 



further discussed in Sect. 3.2. 



The terms involving scalar and pseudoscalar resonance exchange 
will also contain largely unconstrained couplings. At chiral order e^p^, 
there are contributions from mass renormalizations on external legs 
and also from vertex-like corrections. The net effect at this order 
turns out to vanish. 

2. The low-energy parts of meson loop diagrams (Fig. ^(c)): 

We refer to these as the unitarity contributions. They constitute a 
genuine low-Q^ effect distinct from that of the resonance component. 
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In Sect. 3.3 we shall describe such unitarity terms and provide a natural 
extension to all scales, without introducing new parameters. 



Before proceding to a description of the calculation, we introduce a pa- 
rameterization in terms of reduced amplitudes {Cj} and {W^j}, 

with 

= / dQ" WiiQ") (21) 

JO 

and r/_| = t/qo = and 7/+o = 1- In addition, we partition each C 

amplitude into additive components as 

^(eV) ^ ^(mtchg) ^ ^(unty) _ ^^2) 

The matching component ^ encompassing the sum of the Born + form 

factor and resonance contributions, is discussed in Sects. 3.1,3.2 whereas the 
unitarity component c'^^^^^^^ is discussed in Sect. 3.3. The contribution of 
each component to the full amplitude is given in Table 1 (c/Sect. 4). 



3.1 Born and Resonance Diagrams 

The class of diagrams involving exchanges of the ground state pseudoscalar 
mesons and of the low- lying spin-one, spin-zero resonances generates con- 
tributions at order e^]P and at order e^]? . However, we already know the 
e^]P contributions because chiral symmetry relates the K — > tttt amplitudes 
to the K'^ — > vr"*" matrix element and we have calculated this in Ref. [^. 
Therefore we focus on the e^p^ piece in the following. We treat first in some 
detail the Born contributions and their corrections which arise from the in- 
sertion of meson form factors. Then we describe the parameter-dependent 
spin-one resonance terms and finally the spin-zero resonance terms. 

3.1.1 Born and Form Factor Contributions 

The Born diagrams do not contribute to M^oo(Q^) (which involves only 
neutral particles) while giving nonzero contributions to both M^+o(Q^) and 
W+_(Q2). For W+o(Q2) we find 

W^^{Q') = ^J{Q\Ml) , (23) 
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with 



6m^ 



1 + 4 



1 + 6 



.m 



(24) 



This contribution is logarithmicahy divergent at high and has an infrared 
1/Q integrable singularity at = 0. In addition, it is suppressed by a factor 
of M^/M^. This suppression is 'accidental' in that it is not required by any 
symmetry at moderate or high values of Q^. The result is shown as the 
dashed line in Fig. 




Figure 4: Born (dashed) and Born plus form factor (solid) contributions. 

The Born contribution to (Q^) is analytically more involved. Once 

one extracts the infrared divergent singularity j^, it reads 

W+_(Q2) = C{Q^) + SiiQ^M^) - S2{Q\Ml) . (25) 

The functions C^Q"^), Si{Q'^ ,m^) and S2{Q'^ ,m^) are given in Appendix^, 
and we display (Q^) ^ the dashed line in Fig. |5[ Again, this con- 
tribution is logarithmically divergent at high Q^. The cusp is due to the 
singularity related to the coulombic rescattering. 

The set of Born diagrams, required by chiral symmetry, provides a good 
description of the very low region, in which the photon 'sees' only point- 
like pseudoscalars. As increases this is no longer true, and one needs 
to account for structure dependence in the couplings. In our model this is 
accomplished by introducing the low-lying resonances. 

We consider first the diagrams involving pion and kaon electromagnetic 
form factors (saturated in this model by the vector meson resonances) . This 
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Figure 5: Born (dashed) and Born plus form factor (solid) contributions. 

is a subclass of all diagrams required by chiral symmetry but has some nice 
features. It does not introduce any new parameters and improves the high- 
behavior of the {M^j(Q^)} while having minimal effect on the model- 
independent Born contributions at low . The results of this improved 
description are shown graphically in Figs. ^ and |5| (solid lines). The anlytical 
expressions are 



W 



+0 



Ml 



Q2 + M2 



Ml 



J{Q\Ml) + 



g2 + Mj 



J(Q\Ml] 



(26) 



and 



Ml 



Q2 + M2 



c(g') + 



Ml 



Ml 



Q2 + M2 I Q2 + M2 



(27) 

where J{Q'^,m^) is defined in Appendix]^. Note that the new contribution 
to W+o is not suppressed by M^/M^ and thus gives a substantial correc- 
tion to the Born amplitude. In the case of , however, the form factor 

contribution has simply the effect of softening the high-Q^ behavior. 

In principle, given the convergence properties of the 'Born -|- form factor' 
contributions, their integrations over can be performed up to infinity. 
This contribution is dominated by the low and intermediate energy regions, 
where the formalism is valid. This gives a first clean contribution to the 
{Ci} coefficients beyond the Born approximation. 
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3.1.2 Resonance Contributions 



Our analysis contains two classes of resonances contributions, spin-one and 
spin-zero. We consider each one separately in the following. It turns out 
that the spin-zero contributions sum to zero, so that only the spin-one con- 
tributions are subject to the matching procedure of Sect. 3.2. 

As noted earlier, chiral symmetry requires the presence of all possible 
vector and axialvector resonance exchange diagrams. In principle these in- 
troduce to the {Wi{Q'^)} a dependence on eight new parameters, describing 
the weak couplings of resonances. Since the analytical expressions for this 
large class of contributions are rather lengthy and do not illuminate the 
underlying physics, we refrain from reporting them here. The only feature 
relevant for our discussion is the general form, 

W,(Q2) = V gafL'HQ') ■ (28) 



a 



Explicit calculation shows that the physical amplitudes actually depend only 
on the seven parameters, 

.(V) -(A) (A) (A) , . 

5l,2,5,6' 55,6 1 92 - 9l ■ l^yj 

Let us consider the high-Q^ behavior of the functions fa\Q'^) appearing in 
Eq. ( p^ ) . Many of them go to a constant at high or fall off as 1/Q^ , thus 
leading to divergences in the integration process. Such behavior has already 
been observed in similar calculations of the electromagnetic mass shift of the 
kaon [^, 0] . This simply means that the resonance dominance approximation 
can be trusted only up to some intermediate energy region and cannot be 
extended up to Q'^ ~^ oo. In Sect. 3^ we shall try to solve both these 
problems (proliferation of unknown parameters and high-Q'^ divergences) 
by requiring that the resonance amplitude contribution match the high-Q^ 
behavior of the {Wi{Q'^)}. 

We consider next the spin-zero resonance contributions. In the absence of 
electromagnetism, the tree level exchange of the scalar and pseudoscalar res- 
onances contributes a major part of the K — > vrvr amplitudes at chiral order 
p'^. Dressing these diagrams with one virtual photon generates contribu- 
tions to the amplitude 5Af^ at orders e^p^ and e^p^. It is easy to convince 
oneself that diagrams with vertices coming from mass matrix lagrangians, 
having already four powers of the pseudoscalar masses, will contribute at 
order e^p^ to On the other hand, diagrams involving derivative ver- 

tices can give rise to effects of order e^p^, which we are interested in. This 
happens through two classes of contributions: 
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1. mass renormalization on external legs, and 

2. vertex correction diagrams, with virtual photons inserted according to 
minimal coupling. 

Upon explicitly identifying and calculating these diagrams, we find an exact 
cancellation between the two contributions. This is identical in nature to the 
one found in Ref. Q for the Born contributions at order e^p'^ . The explicit 
results (showing the cancellation) can be found in Appendix]^. 

3.2 The Matching Procedure 

As stated in the above discussions, the resonance exhange contribution pro- 
vides a good description for the {W^i(Q^)} only up to some intermediate 

region, beyond which the quark electroweak and strong interactions pro- 
vide the correct framework. Experience in similar hadronic calculations has 
shown that the transition or matching region occurs for Q between 1.5 GeV 
and 2 GeV (or 2 < (GeV^) < 4). The genuine short distance contri- 
butions were studied in the chiral limit in Ref. Q. The outcome was that 
the short distance contribution to the {Wi{Q'^)^ is rather small compared 
to the long distance component. Corrections to the chiral limit cannot dra- 
matically change this qualitative picture. We can imagine assigning a 100% 
uncertainty to the short distance component around the central value given 
by the chiral limit calculation. Even in this case the long distance contribu- 
tion would dominate and our ignorance of short distance physics would not 
significantly alter the final answer. For our purposes, the most important 
feature emerging from this analysis is that for > the {Wi{Q'^)} can 
be set to zero, even if we do not know the details of this transition. 

On the other hand, for low and intermediate we have reliable ex- 
pressions for the {Wi{Q^)}, i.e. the most general parametrization implied 
by chiral symmetry and the low-lying part of the hadronic spectrum. The 
only problem with these expressions is the presence of a large number of 
resonance parameters unconstrained by phenomenology. In what follows 
we shall present a set of reasonable theoretical constraints to be imposed 
on them. The underlying strategy is to use on the one hand the few phe- 
nomenological inputs presently available, and on the other, to enforce the 
transition to the high-Q^ region, meaning in our case that the {Wi{Q'^)} 
have to approach zero in the matching region. 
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3.2.1 Physical Constraints on the g^^'^^ 

Although the {g^^'"^^} of Eq. ( [T^ ) (see also Eq. (p^)) are not predictable 
from a purely theoretical approach, some information can be gleaned from 
the phenomenology of kaon decays and assorted theoretical requirements. 

The phenomenology of kaon decays, especially the radiative kaon decays, 
allows in principle the extraction of a large number of order p'^ constants 
of the weak chiral lagrangian Q. Assuming resonance dominance for these 
couplings (or whenever possible, subtracting the short distance contribution) 
allows one to extract information on the resonance coupling constants. The 
present experimental situation does not, however, yet permit a complete 
implementation of this program, as only limited information is available. 
From K 27r, Svr data and assuming resonance saturation of the relevant 
0{p'^) counterterms, one finds 

45^-5^= 0.43 , (30) 

with a 20% uncertainty associated with the extraction of 0{p'^) coupling 
constants. |5| The — >■ tt transition provides additional information. 
The decay amplitude depends on a parameter w^, whose experimental value 
is w+ = 0.89^014. It receives both long and short distance contributions. 
Using resonance saturation at fichPT = Mp and including explicitly the 
penguin contribution we find 

^'^+ - log - TTTT, - log TT 



V2 



^6 



2-g): >) - FvGvti 



(31) 



At present no other phenomenological constraints are available and we thus 
turn to the description of the theoretical ones. 

In the first place there are two conditions coming from the analysis per- 
formed in the chiral limit |^. Let us recall the reasoning behind this. In the 
chiral limit, using soft-pion methods, one can relate the K — > vrvr amplitude 
to the off-shell K'^-to-'K'^ matrix element. Moreover the invariant amplitude 
Aii+^T^+ is expressible as 

Ak+^^+ = / dQ^ (Q2) . (32) 
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We get one condition by demanding that ^4++ (Q^) vanish at infinity (no 
quadratic divergences). A second condition comes from demanding that 
(Q^) have no short distance component, i.e. that it vanish in the match- 
ing region defined above. Together these amount to 

hm (Q2) = and A++ {^?) = . (33) 

Here we have introduced the matching scale // and according to our previ- 
ous discussion, let it vary between 1.5 GeV and 2 GeV. The following two 
constraints then emerge from Eq. (|33[), 

-(V) , -(V) f ( \ ^ ^'^ 
9i +92 = fviP') 



Fi - Fi 



^i^) - 3#^ = /^(/i) = ^ ^-^0 T . (34) 

^ ^ ■ Ml + ^2 

Adopting the same argument, we require that Wj^ (Q^) and W^oo(Q^) van- 
ish in the matching region, 

W+- (/i^) = and W+Q (/i^) = . (35) 

We do not include the analogous condition for M^oo(Q^) because this func- 
tion, independent of any choice of the parameters, is already very small in 
the matching region (it does not contain any term going to a constant for 
high We believe that the conditions in Eqs. (|3^),(||),(|3D,(|5D form 

a consistent set of physical requirements and provide us with a solid ba- 
sis for any attempt to obtain a sensible answer for the so-called matching 
amplitudes {C J'^^^^'^) } . 



3.2.2 Results 

The above constraints are well motivated and reasonable, but are not suf- 
ficient to completely fix all of the resonance parameters. At this stage, we 
could use specific models of resonance behavior to estimate the remaining 
parameters and then accept the range of model dependence as an estimate of 
our uncertainty. In doing so, however, we have found that models generally 
give a rather small range of results. (The exception concerns 

^(mtchg)^^ The 
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Figure 6: Probability density function for C. 



(mtchg) 



reason is that the matching constraint is more important than the remain- 
ing parameters. Therefore, rather than using particular models we elect 
to follow a more model independent procedure of allowing these remaining 
parameters to vary completely over their reasonable physical range, and to 
use the resulting variation to estimate the error bars for our result. 



The conditions described above imply a set of linear equations for the 

parameters gl^''^\ In particular, we can express all the parameters in terms 
of just two coupling constants. That is, using what we believe are well 
founded physical constraints, we select a two-dimensional hyperplane in the 
parameter space which we call the reduced parameter space. We chose as 
independent variables spanning this plane x = and y = — gi^^ ■ 
By looking at the structure of the constraints one discovers that the other 
parameters depend on x,y and fi as 



(a;,|u). 



and g^'^^{x,y,ii) 



(36) 
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Figure 7: Probability density function for Qq . 

We are now in a position to determine the component C^^^'^^^^ of the 
full amplitude Cj which is determined via matching. The construction de- 
scribed above allows us to express the predictions for each d^^^^^^ as a linear 
function of x and y, 

(12.2 - 0.72x + 0.02y) ± | - 1.0 + 1.35x + 0.006y| , 
(-9.2 + 11.47x + 3.2y) ± | - 4.3 + 6.O2; + 1.0y| , 
(-0.0035 - 1.53x) ± I - 0.0003 - 0.42x| . (37) 

The uncertainties cited in Eq. (^) are associated with the matching proce- 
dure (the variation of the parameter fj.) . Still, this leaves freedom to pick any 
value for (x, y). We can further narrow our predictions by requiring that all 
the couplings simultaneously (as a function of {x,y)) have a 'natural' order 
of magnitude, which can be shown to be 0(1). The existence of a region 
in the (x, y) plane such that this happens is not guaranteed a priori and 
provides a good consistency check for our method. We call this the physical 
region in the reduced parameter space. Studying the explicit dependence of 
the parameters on x, y and /i we are lead to define the physical region as 
X : 0.5 — > 1.5 and y : —1 — > 1. 



C 



(mtchg) 
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Figure 8: Probability density function for c^^^^^ 



A first qualitative conclusion can be already drawn by looking at Eq. (|37| ) 
with (x, y) restricted to the physical region and no further assumptions. The 
expressions show that C^^^^^ depends very weakly on the choice of {x, y) 
in this region, and thus we arrive at a good prediction for this parameter. 
Cq™**^^^^ has a moderate dependence on x and does not depend at all on y. 
This implies that the decay amplitudes can be predicted in our model 
with a reasonably small uncertainty. Problems arise in the expression for 
^(mtchg)^ which displays a fairly strong dependence on x and a moderate 
one on y. In this case, even confining ourselves to the physical region we 
obtain a spread in the answers of about 100%. The only definite prediction 
emerging is that this contribution is not big. 

Quantitative estimates for our results and the attendant uncertainties 
can be obtained by constructing probability distributions for the {Cj-™*'^'^^'*}, 
by means of a survey of the parameter space. We scan the region defined 
by {—3 < X < 3,-3 < y < 3} using gaussian distributions for the input 
parameters. The choice of the parameters in the distributions is made in 
such a way to enhance contributions coming from the physical region. In 
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^ ^ :^ 

ia) (b) (c) 

Figure 9: Two-photon insertions: (a) T^^, (b) V^^, (c) S^^. 

view of this, we choose the central values as x = 0.8 and y = and set 
variances equal to 0.4. The uncertainites cited in our results correspond to 
a 68% probability. Results for the are given in the second row of 

Table 1. 

3.3 The Unitarity Diagrams 

Next we discuss in detail the class of diagrams schematically represented 
in Fig. 3(c). We are interested in the nonlocal part of these diagrams, 
representing the genuine propagation of mesons at low energy. The high 
momentum part of these diagrams produces (on general grounds local 
effects that can be reabsorbed into the definition of the 0{p^) low energy 
constants. In our approach, however, the local component at 0{p'^) is im- 
plicitly contained in the resonance exchange diagrams and would show up 
explicitly upon expanding the resonance propagators. Keeping only the low 
momentum part of the meson loops ensures that the different contributions 
we are including in our calculation do not lead to double counting. Since 
the separation of local and nonlocal components in the meson loop diagrams 
is not free of ambiguity, we shall be careful to describe and motivate our 
prescription in the following. 

3.3.1 Identification of the 0{e'^p'^) Contribution 

Our task in the following is to identify the part in each meson loop diagram 
which, upon contracting the photon legs, will lead to ©(e^p^) contributions. 
The loop contributions to W^i, can be obtained by starting with any meson 
loop diagram which contributes to K ^ vrvr and attaching two photons in 
all possible ways. We focus first on the subclass of diagrams obtained by 
attaching the following two-photon insertions, 

T^,u{p,q), Vf,u{Pi,q), Sf,u{pi,q), (38) 
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as represented in Fig. ^ 

For definiteness let us refer to the bare topology of Fig. In this case 
one can insert T^uip, q) on internal and external legs, Vf^u{pi, q) in the weak 
vertex and S^i,{pi,q) in the strong vertex. The external leg insertions will 
generate wavefunction and mass renormalizations. The other insertions will 
give rise to diagrams like 



Dr. 



D^fjcak 



<fp 



Vw{pi) Vs{pi) 



{k-pf 



Vs{pi 



p^ 



{k -py 



strong 



d^p 



Ml, 



Vw{pi 



{k -pY' 



p^ 



A D^^''{q)T^,{p,q) , 

A D^-'iq) S,.{p,q) . 

(39) 



We wish to isolate the dominant contributions at low momentum (small 
p). Therefore we Taylor expand each tensor insertion around pi = (in 
addition we must expand each coefficient of the Taylor series in powers of the 
pseudoscalar meson masses; for notational convenience we don't explicitly 
display this step). Considering for example the self-energy insertion, one 
has 



dT y 1 

dp a 2! ' dpadpp 



(o,g) + 



(40) 



Evaluation of the integrals in Eq. (|39| ) can be done term-by-term in the 
series. The analysis of each term is very simple. The tensor structure 
d"- /dp"'T^^{p, q)\p=Q factorizes out of the integration over p, and the two- 
photon insertion is replaced by a meson vertex of order p" . This makes 
power counting transparent — after contracting the photon legs it is easy 
to realize that only the first term in the above expansion produces an effect 
of order e^p^ in the kaon ampliude. In the example considered one has 



d^p 



VwiPi) Vs{pi 



{k -p)' 



p^ 



d\D^^%q)T,,{Q,q) +. 



(41) 

This procedure allows us to identify and interpret the relevant contri- 
butions at order e^p^. The integral of T^,y(0, g), weighted by the photon 
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P2(P) 



K(k) 




Figure 10: Loop diagram with internal particles Pi and P2- 



propagator in the above expression, is exactly the expression for the elec- 
tromagnetic self-energy of a charged meson in the chiral limit. Thus the 
insertion of T^jy in a loop diagram reproduces the effect of inserting the 
electromagnetic mass difference into such a loop. Eq. ( |4l| ) then represents 
the meson diagram of Fig. IC with a mass shift insertion on the P2 inter- 
mediate leg. Analogously, insertions of V^^, and S^^, yield Fig. IC but with 
the weak and strong vertices replaced (respectively) by constant vertices of 
order e^p^ , proportional to the chiral couplings Qemw and gems [01 • In other 
words these contributions are the counterparts to what were called implicit 
diagrams in the ChPT calculation of Ref. Q. Their presence in the disper- 
sive matching model is welcome because they provide imaginary parts to 
the amplitudes, ensuring at this order the behavior required by unitarity. 
These expressions when pcirajnietrized in terms of Qemw 

and 5M^, are iden- 
tical to those obtained in ChPT (see Eq. (30) of Ref. Q). Corresponding 
to these contributions will, of course, also be nonvanishing real parts, whose 
treatment is the subject of the next subsection. 

What becomes of the class of diagrams having separated photon vertices? 
The basic result is that they start contributing to K ^ vrvr amplitudes 
at order e^p^. A short argument for this is as follows. Upon contracting 
the photon legs, it is easy to recognize that these loop diagrams have the 
following peculiarity — their intermediate states always involve a photon 
(they contain one photon plus one or two pseudoscalar mesons). Let us now 
consider the diagrams as analytic functions of the external four momenta and 
analyze their imaginary parts as obtained by using the cutting equations. 
The above observation on the structure of the intermediate states, together 
with the form of the lowest order vertices and the phase space, implies that 
the imaginary part of these diagrams is of order e^p^. The long distance 
portion of the real part of the loop, which is all that we are interested in 
here, will then appear at the same chiral order. 
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3.3.2 Real-parts of the Unitary Amplitudes 

In the previous subsection we sliowed tliat tlie relevant part of the unitarity 
contributions at order e^p^ can be calculated with a simple recipe: the 
photon insertion factorizes out and one is left with the calculation of meson 
loop integrals with mass insertions on internal or external legs and weak 
or strong vertices replaced by constant vertices. As we have stated in the 
introduction to this section, we want to keep only the low energy part of these 
meson loops, the one that cannot be mimicked by any local counterterm or 
resonance-exchange diagram. We can best describe this procedure using a 
simple case, the two-pion loop, which also turns out to be the most relevant 
for the physics. The extension to all other diagrams is then straightforward. 

The basic function entering the description of two-pion loops is Jtt-w{s), 
which is given in dimensional regularization by 



where 



fd-i 

L>e + log— J + 1 +/31og ' 



My ' "Kp + i 



(42) 



De^[j^-l + \ogA7i + lj , (43) 

is the scale parameter introduced in dimensional regularization and /3 
is the pion velocity in the kaon rest frame. The divergent piece and the 
scale dependent logarithm in Eq. (^) are clearly local effects. On the other 
hand, the last term and the log term are associated with the low energy 
meson propagation. Finally, an explicit cutoff calculation shows that the 
additive factor of one has to be included in the long-distance part. These 
considerations lead us to introduce a separation scale such the short and 
long distance parts are defined as 

Jif (a4) = ^ (d,+ log 4 I ■ (44) 



There is an inherent ambiguity in the separation scale which cannot 
realistically be assigned a unique value. Therefore we let it range between 
Mk and Mp, associating the corresponding variation in the result as the 
theoretical uncertainty. These unitarity corrections come with moderate 
uncertainties except for the case of Cqo (c/ Table 1). 
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3.4 EM Corrections to the Isospin Amplitudes 

Let us consider some phenomenological consequences of our analysis. We 
refer the reader to Sect. 2 and to Sect. 4.3 of Ref. for an introduction to 
formalism used in the following. In the presence of electromagnetism, the 
amplitudes involving charged particles (^-| and ^+o) contain infrared sin- 
gularities. For each such amplitude, the infrared singularity can, on general 
grounds, be isolated in an exponential factor that multiplies an infrared- 
finite amplitude which can itself be expressed as an expansion in powers 
of alpha. Upon considering the emission of soft photons with energy up 
to some experimental scale iv, the infrared divergences disappear from the 

decay rate expressions, leaving w-dependent factors G-i (cj), G+o{uj) which 

multiply the square moduli of the infrared- finite amplitudes. This process 
has been explicitly described in Sect. 4.3 in Ref. Q for the K'^ — > vr+vr^ 
mode. 

Starting from the infrared-finite amplitudes in the charge basis, we can 
define the would-be isospin amplitudes from the following linear combina- 
tions, 

2 . 1 . 
Ao = 3-^+- + 3"^oo ' 

A2 = ^M+--Ao) , (46) 
A2 = 3-^+0 • 

In the absence of electromagnetism and any other isospin breaking interac- 
tion, we then have A2 = , and the amplitudes of Eq. ( ^6[ ) truly represent 
transitions to pure isospin states. Using the same logic one can perform an 
analysis of the unitarity condition ^ , leading to the following parametriza- 
tion of the K vrvr infrared finite amplitudes, 

v2 

Aoo = (^o + <5^r)e^(''«+^o)-\/2(^2 + Mr)e'^'^'^^'^ > (47) 



2 

The calculation performed in this paper gives us knowledge of the 
We find the shifts in the isospin amplidudes to be 



Mr = 4^ + 3C00 ) = (0.0253 ± 0.0072) IQ-'M^^o 
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= ^F^lf ^ ^ - Coo) = (0.0118 ± 0.0063)10-^ M^o , 
M^"'" = p^^^C+o = -(0.0080 ±0.0088)10-^ M^o . (48) 

In our numerical evaluation we have used a value for obtained from a fit 
to data not including radiative corrections. |^] This introduces an ambiguity 
in of order a which affects at order and thus is beyond the 

accuracy we are working at. As a byproduct we obtain also the effective 
A/ = 5/2 amphtude, 

^5/2 = — (C+- - '^00 - C+O) 

= (0.0137 ± 0.0097) • 10^'^M^o . (49) 



4 Conclusions 

The problem of determining electromagnetic corrections to nonleptonic kaon 
decay is a formidable one and has long resisted understanding. In this paper 
we have employed a 'dispersive matching' approach which provides a frame- 
work that is, in principle, general and model- independent. This dispersive 
setting was first advocated by Cottingham ||lO| and has been recently em- 
ployed in Ref. At a practical level, however, a rigorous implementation of 
this program is plagued by a lack of sufficient input data. We have been able 
to overcome this obstacle by pointing out (on rather general grounds) how 
long range and intermediate range processes are expected to dominate the 
physics and then performing a calculation which incorporates the relevant 
ground state and resonance degrees of freedom. All possible tree-level am- 
plitudes and a subset of loop amplitudes are taken into account. The latter 
component ensures that our amplitudes have the imaginary parts required 
by unitarity. 



Table 1: The {Cj} Amplitudes 





C+- 


Coo 


e 


fO 




e p' 


-3.3 ± 1.7 





-3.3 


± 


1.7 


e^p^ (Matching) 


11.6 ±0.3 


-1.3 ±0.7 


-0.4 


± 


4.5 


e^p^ (Unitarity) 


6.5 ± 1.5 


3.1 ± 1.4 


-3.4 


± 


1.2 


Total 


14.8 ±3.5 


1.8 ± 2.1 


-7.1 


± 


7.4 


ChPT 


14.1 ± 12.5 


0.9 ±6.7 


-4.2 


± 


4.6 
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Contributions to the {Cj} are shown in Table 1. The first row displays 

terms of order e^p^ calculated in Ref. . The equal values for and C+o 

are due to the absence of a A/ = 5/2 component at lowest order. The second 
and third rows display terms at order e^p^ , arising from the analyses done in 
Sects. |3.1|-3.2 and in Sect. 51^. The fourth row summarizes the total result 



obtained within the dispersive matching approach. For comparisons sake, 
we also cite the ChPT results in the final row. We obtain EM corrections 
whose central values are in reasonable accord with our earlier ChPT calcu- 
lation [Q] but whose theoretical uncertainties are substantially smaller for 
the vr+vr^ and ir^ir^ modes. Only the — > Tr+Tr*^ determina- 

tion produces a less precise value. It is not hard to recognize the reason for 

this. In our language, C-| is dominated by the low-Q^ Born contribution. 

The resonance contribution, important at intermediate Q^, introduces only 
a moderate uncertainty. On the other hand, for C+o the low-Q^ contribu- 
tion is very small, being suppressed by a factor M^/M^ (see Sect. 3T). C+o 
is thus dominated by intermediate effects, which are plagued by a sub- 
stantial uncertainty that our constraints have not completely eliminated. 
We could turn to model-dependent frameworks to attempt to narrow the 
quoted error bars. However, this apparent improvement would likely be il- 
lusory, since our understanding of models is too weak for any specific model 
to be trusted in a calcuation such as this. Thus we feel that our quoted 
error bars are a reasonable measure of present uncertainties. Note however, 
that the uncertainty in C+o is not much of a problem because of the overall 
smallness of the effect. 

A key result of our calculations is that the electromagnetic corrections 
to the weak amplitudes are smaller than naive estimates might indicate. 
Part of the reason is the partial cancellation in the leading chiral transition 
that we detailed in Ref. ||^]. In addition, only about a third of the overall 
electromagnetic effect goes into a modification of the 1=2 final state - the rest 
is harmless as it contributes to the much larger 1=0 final state amplitude. 
Although the work done here constitutes a crucial step in our study of EM 
corrections to nonleptonic kaon decay, there remain several additional issues 
which we shall address in a future publication. |^] Chief among these is how 
to correctly extract the electromagnetically corrected K — > vrvr amplitudes 
from experimental data. We shall discuss the underlying theory in some 
detail, as well as suggesting the proper procedure to be followed in the 
experimental analysis. Another topic to be covered, of great current interest 
in studies of CP violation, involves the ratio e'/e. The calculation done here 
leads to a value for the EM correction to e'/e (commonly denoted as r^EM)- 
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We shall also provide an improved determination of the phase of e'/e. 

The research described here was supported in part by the National Sci- 
ence Foundation. One of us (V.C.) acknowledges support from M.U.R.S.T. 
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A Loop functions 

In this appendix we give the analytic form of the functions entering in the 

Born term of (Q^)- It is convenient to express them in terms of four 

simpler functions arising from the integration over the angular variables, 



The Fi are 



= Q' 



Q^F4(Q^)+Fl(Q^Mi) 
2Ml 



2(Fl(Q^m2)-/32F2(Q^m2)) , 
2p^[Fi{Q^,m^)-F2{Q\m^)) , 
2Fl(Q^m2)+F2(Q^m2) . 



(50) 



1 



2m? 
8m4 



1 + 



2m^ 



1-4 



1 



2M|g2 I 



log 



1-/? 

1 + 



log 



1-13x1 + 



g2 + 4M2 



(51) 



B Scalar and pseudoscalar resonance contribution 

Let us define the parameter /em as follows, 



-^EM = -T" / dQ 
47r Jo 



(52) 



The scalar and pseudoscalar resonance contribution to the vertex-correction 
diagrams at order e^p^ is given by the following expressions for A+-, 



A 



■+- 



Iem 



c., (2a^ I " ^ 
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V2 



Cd 



ml - Ml)gf + 3(Mi + -Ml)gl 



Ml 



F^Fk 



92 Cr, 



4M|- + 2M2 



Ml 



+ 



4^2 .s~ 
9iCd 



Ml 



Ml 



F^F 



K 



and for A+o, 



Ml 



(53) 



3F^F, 



K 



Cm 54 



sMl 



Ml 



Ml 



+ 



F^Fk 
1 

'FfFl 



Cd 



Mhf + {Ml + -Ml)gl + Mlgl 



92 Cm 



2Ml 



Ml 



Ml 



4 Ml 

F^Fk M|, 



(54) 



The mass renormalization effect is given by the negative of this expression, 
with Jem replaced by ^M^. However, the long distance contribution of SM"^ 
is just given by /^m and thus these terms cancel each other. We neglect any 
residual intermediate energy component that may occur. 
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